Ionization potentials in the limit of large atomic number 
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By extrapolating the energies of non-relativistic atoms and their ions with up to 3000 electrons 
within Kohn-Sham density functional theory, we find that the ionization potential remains finite 
and increases across a row, even as Z — >■ oo. The local density approximation becomes chemically 
accurate (and possibly exact) in some cases. Extended Thomas-Fermi theory matches the shell- 
average of both the ionization potential and density change. Exact results are given in the limit of 
weak electron-electron repulsion. 
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A central problem of electronic structure is the calcu- 
lation of the ground-state energy of the electrons of any 
atom, molecule, or solid, within the non-relativistic Born- 
Oppenheimer limit. Density functional theory (DFT) is 
a popular choice, balancing computational efficiency with 
useful accuracy. The original DFT was that of Thomas 
[H and Fermi j2] , TF theory, in which a local density ap- 
proximation is made for the the kinetic energy and the 
electron-electron repulsion is approximated by the sim- 
ple Coulomb energy of the charge density. In the 1970's, 
Lieb and co-workers Q showed that the TF energy be- 
comes relatively exact for neutral matter as Z — > cx) in a 
specific way. The energy, E, grows in magnitude as Z'^^^, 
where Z is the total charge. For atoms and their ions, 
the leading corrections in powers of Z^^^^ were found 
by Scott Q, Dirac @, Schwinger and others d, 0], as 
summarized by Englert Q. These corrections are given 
exactly by extended Thomas-Fermi (ETF) theory, which 
includes both the gradient correction for the kinetic en- 
ergy (one-ninth the von Weisacker functional Q), and the 
local density approximation for exchange (LSDX [ic!]). 

However, TF theory and its extensions are insuffi- 
ciently accurate to predict chemical properties jlljl . Mod- 
ern DFT uses the Kohn-Sham (KS) scheme, in which only 
a very small fraction of the total energy, the exchange- 
correlation (XC), needs be approximated. But the idea 
of asymptotic correctness was recently extended to KS, 
relating the success of exchange GGA's such as PBE fl2l| 
for total energies to their recovery of the (Z^ ) term in 

[li]. The re- 
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the expansion of the exchange energy 
lation between semiclassical and local density approxi- 
mations ^5i| contributed to the creation of PBEsol 16], 
a PBE-like functional that is nearly optimum for solids 
near equilibrium but not for atoms and molecules, and 
to revTPSS 17], a nearly optimal semi-local functional 
for all three kinds of systems. 

But total electronic energies are irrelevant to chem- 
istry. Only differences matter, such as the ionization po- 
tential of an atom (/ is the energy difference between 
the positive ion and the neutral) or the dissociation en- 
ergy of a chemical bond. How relevant are asymptotic 
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groups of the periodic table (from all<ali metals to noble atoms) 

FIG. 1: Ionization potentials of the main groups in the limit of 
large row number of the periodic table, calculated using exact 
exchange, the local (spin) density approximation (LSD), and 
PBE; ETF denotes extended Thomas-Fermi theory. 



expansions for these quantities? The asymptotic expan- 
sion for E is in powers of Z~^l'^ , so if / remains finite as 
Z ^ CO, the neutral and ion energies must agree for the 
first seven powers in such an expansion, a truly remark- 
able balancing act between quantum effects, the Pauli 
principle, and the Coulomb forces of nuclear attraction 
and inter-electron repulsion. In this letter, we demon- 
strate by both calculation and analysis that (i) / has 
no limit as Z — ?> oo, but remains column-dependent (ii) 
that each column has a finite limit; (iii) the local (spin) 
density approximation (LSD [10]) of KS theory becomes 
very accurate (if not exact) for / for certain cases; (iv) 
ETF theory becomes very accurate (if not exact) for the 
average of / over an entire shell; (v) the shell-averaged 
difference in density between the neutral and its ion ap- 
proches that of TF. We demonstrate these statements in 
the limit of weak interelectron repulsion. 

Our most important results are shown in Fig. [TJ We 
plot / from various calculations, extrapolated to infinite 
row number, versus the column number for main-group 
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elements (s and p valence shells). We calculate exchange 
exactly [3 using the optimized effective poten- 

tial (OEP, which here should be indistinguishable from 
Hartree-Fock '2^), extrapolating all values to Z — > 00. 
At the exchange level, LSD and PBE are almost exact 
for p- valence elements, and are highly accurate but in- 
exact for the s-valence cases. Furthermore, ETF yields 
Q a single number (3.15 eV), very close to the s- and p- 
average (3.02 eV). When correlation is included, gradient 
effects are slight, and it is in the regime of large electron 
number that approximate density functionals work best, 
sometimes exactly 13, IS IHj ■ We speculate that LSDX 
on accurate densities becomes almost exact in this limit 
for p-shell cases, that ETF is exact for some shell-average, 
and that our XC results arc extremely accurate and prac- 
tically impossible to calculate with any other method. 

To understand why local functionals become accurate 
in this limit, begin with total energies of neutral atoms, 
whose large- Z expansion is 



(1) 



where Eq{Z) is the energy of an atom with atomic num- 
ber Z and charge Q, and the c^^^ are coefficients depend- 
ing on the degree of ionization, q ~ Q/Z. We use atomic 
units throughout. The neutral coefficients were derived 
via semiclassical analysis by Schwinger and Englert 
The TF energy is exactly — c^°'' Z'^/^. The second term 
[3| comes from the s-orbitals at the nucleus and must be 
treated quantum mechanically. The third term is deriv- 
able in ETF theory of which 2/11 arises from the 
gradient correction to the kinetic energy, and 9/11 from 
LSDX. When treated in terms of the potential, the di- 
vergence at the nucleus is avoided Q. 

The extension of these ideas to / has proven more dif- 
ficult. Terms of higher order than those shown in Eq. (1) 
oscillate with Z, as a precursor to the periodic vari- 
ation of chemical properties that is missed by ETF, but 
well-described in KS DFT. The oscillations in / domi- 
nate over trends with Z~^^^. While numerous studies 
exist in the literature 221 for fixed electron number N 
with Z — > 00, we are interested in I{Z) = Ei(Z) — Eo(Z) 
as Z ^ 00. Within TF theory, Lieb proved [2j| that / 



does not grow with Z, and by considering c^*^^ as g ^ 0, 
Englert showed I^^ SA-^/s/Jq ^ 1.29 eV, where A = 
32.729416 is a known constant [1], and a = (97r2/128)i/3. 
Even this simple result requires explanation, because /i, 
the chemical potential, is zero for the neutral atom in 
TF theory, suggesting / should be too. But the TF en- 
ergy is the smooth envelope of Eq{Z) as a function of 9, 
whereas the true energy consists of line segments between 
integer values 2J|. Thus /i = dE/dq — —I for the exact 
system, but the TF energy behaves as g^/"^ for small q. 
So fiTF = 0, but the better value of I^^ is the energy 
difference Q with Q = 1. 

Next we discuss KS DFT, in which the (non- 
interacting) kinetic energy is not approximated, but is 
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FIG. 2: OEP ionization potential / (in eV) versus Z~^^^ for 
main groups of the periodic table. Also shown with green 
lines is the noble atoms LSDX curve. 



found exactly from the KS orbitals. We perform KS 
self-consistent calculations for atoms and ions up to 2938 
electrons using LSD and PBE XC functional approxima- 
tions, as well as the exact OEP exchange. These were 
done using the Engel code fio^, but with tightened con- 
vergence criteria and maximum numbers of orbitals, and 
a logarithmic radial grid with 800 points. In Fig. [2] we 
show / versus Z~^l'^ for each main-group column of the 
periodic table. In all cases, the behavior is almost lin- 
ear as a function of Z^^l'^ for all Z > 169, so we ex- 
trapolated these curves using a parabolic fit in Z~^l^ 
and found the ionization energy for Z — > 00 as shown 
in Fig. [1] The spherical approximations of the density 
(LSD, PBE) and of the potential (OEP), used in the En- 
gel code (see Ref. [l^), give errors less that 0.1 eV for 
/. We use electronic configurations based on the aufbau 
principle and Madelung rule [25|. For the noble gases, 
Z = n[ri} -I- 6n -I- 14)/6 - A{n){n/2 + 1), where n is the 
row number and A(n) = for even and 1 for odd rows. 

To understand in detail the results shown in Fig. [1] 
which are also tabulated in Table I, we begin at the ex- 
change level. Both PBE and LSD exchange are almost 
identical to the OEP values for the p-group elements, 
with a maximum difference between them of 0.02 eV, and 
of cither from OEP of 0.08 cV. This is not so for the alka- 
lis and alkali earths, presumably because they have only 
one or two electrons outside a closed shell, with accom- 
panying self-interaction error of approximate functionals. 
The ionization of p-elements involves removing electrons 
from a full (or almost full) shell with ~ (n -I- 2)^/2 elec- 
trons, where n is the row number for even n, and rows n 
and n + 1 have the same structure. 

In fact, Englert also showed that the TF result is not 
correct as Z 00. The terms of 0{Z'^^'^) in Eq. (1) 
also yield a finite contribution, which is included in ETF, 
yielding an / of 3.15 eV, very close to the average over 
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TABLE I: Extrapolated ionization potentials / (eV) of main 
group elements. Mean absolute difTerences (m.a.d.) are taken 
relative to OEP for X, and PBE for XC. The last two columns 
show the electron affinity A (eV) (estimated as / — l/{r) in 
atomic units) and the average radius (r) (bohr) of the ion- 
ization density, in the Z ^ oo limit, using PBE. For ETF, 
1=3.15 eV, (r) = 5.6 A, and A = 0.58 eV. 




FIG. 3: Exchange contribution, averaged over shell, to ioniza- 
tion potential for Bohr atom with many electrons; blue circles 
are exact, open circles are LSDX on exact density, and black 
dashed line is LSDX on TF density. Solid lines are cubic fits 
to the last ten circles. 



both s and p-shell values (3.02 eV). 

To check this is no accident, consider the simpler sys- 
tem of atoms with an infinitesimal electron-electron re- 
pulsion, A, sometimes called Bohr atoms. The orbitals 
are hydrogenic, requiring no self-consistency and simpli- 
fying the integrals 26[. One finds that l"^^ is exact for 
large Z at A = 0. In Fig. |31 we show the exchange cor- 
rection (divided by A) to / for LSDX applied to the TF 
density (yielding 8{2/3)^^^ / {Stt^) « 0.2360), to the exact 
densities (each averaged over entire shells), and exactly. 
All three match as Z — oo, but a small error remains 
if, e.g., just the s-shell is used. Thus we speculate that, 
for real atoms, LSDX (in a KS calculation) matches the 
average over the entire shell as Z — >■ oo. 
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FIG. 4: q = 1 ionization density, 4ivr^{natom.{r) — nion{r)), as 
Z ^ OQ for the average over the groups III-VIII, for the alkali 
series, and within TF theory. We use LSDX KS densities. 



Next, we discuss the DFT calculations with correla- 
tion, which remains finite as Z ^ oo and varies across 
a row. The differences between PBE and LSD are rela- 
tively small, giving greater confidence in both. The max- 
imum deviation between them for p-elements is 0.14 eV, 
comparable to the deviations of these functionals at the 
exchange-only level from OEP for the alkali and alkali 
earths. Thus the gradient corrections are not vanishing, 
suggesting that while both calculations are accurate, nei- 
ther is exact. The PBE average, 3.61 eV, is our best es- 
timate of a universal ionization potential, defined as the 
limit of / averaged over the n-th shell, as n — > oo. 

The other major descriptor of chemistry is the electron 
affinity A{Z) = Eo{Z) - E^i{Z). Within LSD or PBE, 
the first negative atomic ion of energy E-i{Z) has no 
stable solution, but A{Z) can still be estimated ^2?'] via 
a charged conductor model, in which I — A = \/ {r), and 
(r) is the centroid of the added charge. Define the radial 
ionization density as 



Anji{Z, r) = 47rr^ {no{Z, r) - ni{Z, r)) , 



(2) 



which integrates to 1. Then choose (r) — 
dr rAnR{Z,r). Table I shows PBE Z oo lim- 
its for /, (r), and A. Averaging over s and p, our best 
estimate for a universal value of A is 0.78 eV. 
We next extrapolate the ionization density via 

AnniZ, r) « /3Anfl(Zo, /3r) + jd[AnR{Zn, r)]/dr, (3) 

which correctly integrates over r to 1. Here Zq = 2935, 
/3 = 1 + &(Z-i/3 _ z-i/3), and 7 = c(Z-i/3 _ z"^/'), 
with fit parameters 6 = 5, and c = —2. Finally, we 
also averaged over the 6 p-shell curves, to find the results 
shown in Fig. S) The TF solution for the infinitesimally 
charged ion has a finite size [1, [2^, i.e., 

rc = lim roiZ) = a A^/^ « 9.0588 bohr « 4.8 A. (4) 

Z—>-oo 
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FIG. 5: Same as Fig. [H but for the noble-gas column of the 
periodic table at various finite Z and in the limit Z — >■ oo. 



in that case, the last shell is spread throughout the entire 
atom, and average gradients contributing to ionization 
vanish as Z — )■ oo, which is not true for real atoms. 

Thus we have established that, in the large-Z limit, 
the periodic table becomes 'perfectly periodic. Moreover, 
local approximations appear to become exact, even for 
energy differences that are relatively vanishingly small in 
this limit. These are new, numerically relevant, exact 
conditions that approximate functionals should satisfy. 

All conclusions are based upon numerical calculations 
and extrapolation. Proving them rigorously is a chal- 
lenge to mathematical physics. This work was supported 
by DOE grant DE-FG02-08ER46496 at Irvine, and NSF 
(Grants DMR-0501588 and DMR-0854769) at Tulane. 



Beyond this radius, An^^(r) is just the radial density 
of the neutral, which has reached its asymptotic form, 
decaying as l/r^. The maximum of this curve is about 
0.1830 at r = 8.855 bohr. The agreement between the 
extrapolated p-shell densities and the TF theory is re- 
markably good, but not exact, while the extrapolated 
alkali ionization density is very different. We speculate 
that averaging over an entire shell would yield perfect 
agreement, as we find numerically for the Bohr atom. 

Finally, we justify why such large atomic numbers 
(larger by a factor of 10 than those of Ref. [ll]) 
needed to get these results. Because of the scaling with 
even Z = 125 only makes Z'^l'^ = 0.2, while 
Z > 1000 brings Z^^l"^ below 0.1, making the extrapo- 
lation much more reliable. In Fig. [5l we show accurate 
ionization densities for the 8-th column of the extended 
table at finite Z. The scaling of the TF ionization den- 
sity is quite different from that of the exact solutions: 
Before extrapolation, even at Z = 2935, the TF ioniza- 
tion density agrees much better with that of the alkalis, 
not the p-shell average. For the same reasons, having 
HF energies for only Z < 100, Englert erroneously con- 
cluded that /ETF .^g^g limit of the alkalis, not the 
shell-average (see Fig. 4-8 and its discussion of Ref. Q). 

Thomas-Fermi theory produces the first term of Eq. 
(1) and extended TF yields an average Z — oo limit 
for the ionization energy, but no pejriodic variation of 
chemical properties and no binding [llj] of atoms to form 
molecules or solids. Within non-relativistic KS theory, 
any reasonable approximation to the XC energy with the 
correct uniform-density limit for exchange will produce 
the total-energy expansion of Eq. (1) and a finite column- 
dependent Z — >■ oo limit for the ionization energy. It 
appears that LSD is extremely accurate and possibly ex- 
act in certain cases for /. We have shown this for shell- 
averages in the limit of weak inter-electron repulsion. But 
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